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Abstract. Albert-Laszlo Barabasi introduced a model [ '] which exhibits the bursty 
nature of the arrival times of events in systems determined by decisions of some humans. 
In Barabasi's model tasks are selected to execution according to some rules which depends 
on the priorities of the tasks. In this paper we generalize the selection rule of the A.-L. 
Barabasi priority queuing model. We show that the bursty nature of human behavior 
can be explained by a model where tasks are selected proportional to their priorities. In 
addition, we extend some of Vazquez's heuristic arguments [j. 4] to analytic proofs. 

1. Introduction 

Nowadays, many human-driven phenomena have become momentous in different areas 
of hfe, for example in economics and in social sciences. Classical models of human dynam- 
ics [11, (>] are based on Poisson processes, so in these models the elapsed time between 
two consecutive events, the so-called inter-event time is exponentially distributed. This is 
so because in these models there are many people who has effect on the system. On the 
contrary, when a single person's actions are considered, in many real cases the inter-event 
time has a heavy-tailed distribution [2]. That is the probability that the inter-event time 
is equal k tends to zero with the speed of k~'^ for a 7 > 1. (For precise definition see [7, 
Ch. 1]) These heavy tailed distributions give relatively large weight to long inter-event 
times, and there are bursts of high activity between these inactive intervals. This is in 
sharp contrast with the exponential distribution, which gives low weight to the longer 
inter-event times. 

Albert-Laszlo Barabasi introduced a model [ ] which exhibited the bursty nature of the 
arrival times of events in systems determined by decisions of some humans. He analyzed 
it using simulations and found that it is in good agreement with empirical experiences. 
That is, for example the timing of e-mails sent by a user can be described with Barabasi's 
model: The distribution of the inter-event time between two emails sent by a selected 
user can be approximated by power-law distributions with exponent close to one, which 
is also obtained by the model in the p — )■ 1 limit. Vazquez [I '] was the first one who 
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studied the model analytically. He confirmed Barabasi's results by determining the exact 
distribution of the priorities and the waiting time at stationarity. He found that in the 
p — )■ 1 limit the distribution of the waiting time is close to a power-law distribution with 
exponent one. 

In this paper, motivated by the work of Barabasi ['^] and Vazquez [11], we generalize 
this model to a more general one, where burstiness occurs as well. This generalization 
shows that this power-law decay with exponent 1 can be obtained in sequences of models 
where tasks with higher priority are more and more likely to be chosen. In addition, we 
extend some of Vazquez's heuristic arguments to analytic proofs. 

Barabasi's model is as follows: Somebody has a todo list which always consists of 
exactly L items, each of which has a non-negative priority. In every discrete time step 
(1,2,...) two things happen: 

• A task is selected for execution according to a selection protocol described later 
and leaves the system, 

• another task arrives to the list. 

The priority of all arriving task including the tasks at the 0-th step are i.i.d. random 
variables with an absolutely continuous distribution function R{x) on [0, oo). Let us call 
a task new at time t if it has just arrived to the list. We denote by Nt the priority of the 
new task and let us introduce 

R{x) := ¥{Nt < x). 
The corresponding density function (DF) is 

r{x) := R'{x). 

Throughout the paper we assume that 

r G ^2(0,1). 

In addition, let us call a task old at time t, if it was on the list at time t — 1 at it is 
still there at time t. We denote its priority by Ot and denote its cumulative distribution 
function by 

Ri{x,t) :=P(Ot < x), 

and its density 

^^1(2;,^) := -^Ri{x,t). 

The selection protocol in Barabasi's model is as follows: Independently in each step 
independently, we toss a coin with probability of heads < p < 1. If heads, the task with 
the highest priority is selected for execution and leaves the system. If tails, a task selected 
uniformly from the list leaves the system. In particular, if p = 1, then we use the highest 
priority first selection protocol, and if p = 0, then the selection is completely uniform. Let 
us denote by r the waiting time of a task to for execution, that is, the number of steps 
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between its arrival at and its departure from the system. Using computer simulations, 
Barabasi found that if p — )■ 1 then the waiting time r has a power law tail with exponent 
1 . These asymptotics hold not only for long lists but even if the list consists of two items, 
that is for L = 2. In this case Vazquez [14] proved analytically that 



In [14] the expected value of the waiting time was also investigated. For the case L = 2 
Vazquez proved analytically that 



in stationarity. If p = 1 this means that there is a task stuck in the queue with priority 
0. For the case L > 2 Vazquez conjectured that 

Conjecture 1.1 (Vazquez). 



in stationarity. If p = 1 stationarity means that there are L — 1 tasks stuck. 
Our work is divided into two main parts. 

1.1. Barabasi model. We extend some of the heuristic arguments of [11] related to the 
Barabasi model to analytic proofs in Section 3. We distinguish two cases according to 
p < 1 OT p = 1. 

1.1.1. Barabasi model with p < 1. First we assume that < p < 1. We extend Vazquez's 
completely correct heuristic argument to an analytic proof in Section 3.1. 

1.1.2. Barabasi model with p = 1. Then in Section 3.2 we consider the case when at every 
time step the task with the highest priority is selected. Although the relevant part of 
Vazquez's Conjecture 1.1 holds but the stationary case will not be achieved, since there 
will never be any tasks with priority almost surely. It is of crucial importance to observe 
that in this case at time t the L items on the list are as follows: Beside the newly added 
item we have the L — 1 lowest priority items among all items which have ever added to 
the list. Therefore we can use the theory of the process of records. We study this case for 
L = 2 in Section 3.2, where all tasks remain on the list for one time step except for the 
records. An item is called lower record (from which we usually omit to say the adjective 
lower) if it has the lowest priority at the time when it is added to the list. When such 
a record arrives, then it remains on the list until a new record arrives. The time elapsed 
between the n-th and n + 1-th records called the n-th inter-record time. It is well known 
that the expected value of the inter-record times grow exponentially with n. Therefore 



(1.1) 
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simulations may indicate the false conclusion that one task remains indefinitely in the 
queue [14] [2]. 

1.2. Generalization of the Barabasi model. In Section 4 We study the Barabasi- 
model for general priority-based selection protocols in the case L = 2. We analyze the 
two main characteristics of this model at stationarity, namely the distribution of the 
priorities of tasks on the list and that of r, the waiting time a task spends on the list 
before execution. To explain these with more detail, we need to introduce some new 
notation. At time t + 1 there are two items on the list. The one which has just arrived 
to the list we call new, and the one which was already on the list at time t we call it old 
task. The distribution of the new task is given by R{x). Our goal is to understand the 
distribution of the old task. The selection protocol is described by the function v in the 
following way: 



where v{.,y) is increasing and Vx v{x, .) is decreasing. Furthermore we assume that 
Vx v{x,y) < Ci < 1 with some constant Ci. 

We use Vazquez's notation, and also refer to the results in [14] for the Barabasi case, 
when v{x,y) = ptx>y + (the model introduced by Barabasi [2]). An other natural 
example is as follows. 

Example 1.2. Let v be v{x, y) = + < p < 1 and the priority of the new task 
is uniformly distributed on [c,l], < c < 1. That is, first we toss a biased coin (heads 
with probability p). If heads, we select a task for execution proportionally to the priorities, 
if tails, the selection is done uniformly. This model takes into account not only the order 
of priorities, but also the proportion of those. 

If tasks are interpreted as being competitors and the priorities as levels of talent of these 
competitors, then Example 1.2 can be explained as follows: In a game, every competitor 
entering the game have to play until he wins for the first time. We also know that the 
outcome of the games depends on the level of talents as it is determined by the selection 
protocol of Example 1.2. That is, with probability p the chances to win are proportional 
to the level of talents of the two competitors, and with probability 1 — p the competitors 
win with probability. The a priori distribution of the talent of a player is uniform on 
[c, 1], and the game has been going on for a long time. We are interested in the distribution 
of the level of talent of a player about whom we only know that he has just lost a game. 
Then the corresponding density function ri{x) is illustrated in Section 4.2, see Fig. 2. 

Very useful notations are the following ones. The probability that the new task is 
selected given the old task has priority s is 



v{x,y) := P(new is chosen | new has priority x, old has priority y). 



(1.2) 
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and the probability that the new task is selected given the old task has priority s is 

= / {l-v{s,y))dR,{y,t). 



JO 

With these notations one can easily describe the evolution of R{x,t). Assume that at 
time t + 1 the old task is T*. Then either 

(a) T* was the old task at time t, or 

(b) T* was the new task at time t. 

In the first case (a), the new task had to be selected for execution in step t, so that T* 
remained in the system, and in the other case (b) the old task had to be selected so that 
T* remained in the system. Hence law of total probability gives us 

t + l)= / ri{s,t)q{s)ds + I r{s)qi{s,t)ds . 
Jo Jo 

V V 

case{a) case(b) 

Assuming the system is in stationary state, the probability that the old task is selected 
given the new task has priority s is 

(1.3) gi(s)= I il-v{s,y))dRiiy), 



JO 

and for the stationary CDF Ri{x) and DF ri{x) of the priority of the old task we obtain 
the stationary equation 

px i-x 

(1.4) Ri{x) = / ri{s)q{s)ds+ / r{s)qi{s)ds. 

Jo Jo 

In the very special Barabdsi case Ri{x) can be explicitly computed [ I 1]. However, in the 
general case the function Ri{x) cannot be expressed by a closed formula. We will use 
Hilbert- Schmidt operator techniques to find an approximation of Ri{x). 
As far as the waiting time is concerned, the distribution of r is: 

{I - qi{x))dR{x) ifA; = l, 

Jq - q{x))q{xY~'^dR{x) if A; > 1. 

If /c = 1, then the task has to be executed immediately, so the probability of this is exactly 
1 — ^1(0;) if the task has priority x. Otherwise, if A; > 1, then the event that a task with 
priority x stays exactly k steps on the waiting list is the intersection of k conditionally 
independent events. When the task is added to the list, it must not be executed, this 
happens with probability gi(x), then it becomes old and has to stay k — 2 more steps on 
the list, the probability of this event is g(a;)^~^, and finally, the task has to be executed, 
which event has probability 1 — q{x). 

In the Barabdsi case these integrals can be explicitly computed [1 1], hence we obtain 



;i.5) P(r = k) 



{T = k) 



^(m"'-(¥)"0^ if^>l. 
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In the limit p — 1 it leads to (1.1). 

In the general case only the expected waiting time can be explicitly computed, see 
Section 3.1. 

2. Probabilistic interpretation of Ri{x) in the Barabasi case 

Consider the Barabasi case withp < 1. As we mentioned earlier, Vazquez [11] computed 
Ri{x) in the Barabasi case. In order to give a better understanding of the nature of the 
model, here we provide an alternative, probabilistic way to compute the CDF Ri{x). Part 
of this method will be used when we prove Conjecture 1.1. 

To get the distribution function Ri{x) of the old task in a stationary system, first 
extend the dynamics backwards in time such that the values of the system at —1, —2, . . . 
are defined to maintain stationarity. Further, let us suppose that the process is stationary 
at time T. Starting from this T, we are looking backwards in time and count how many 
other tasks had an old task to compete with in order to be able to stay in the list until 
T. We define the following three events: 

• -Knew := {The selection protocol turns out to be the random choice and 

a new task is executed. } 

In this case the new task leaves the list immediately. Thus, the old task does not 
have to compete with its priority at all, so we do not have to count these cases. 
This event happens in each step with probability 

• -Roid '■= {Random protocol is chosen, and the old task is selected.} 

In this case the arriving new task becomes old and the old task before it leaves the 
system. Since the priority distribution of the new task is just R{x), considering 
the old task's priority distribution, the system restarts at each event of this type. 
This event happens in each step with probability 

• -Rcomp := {Priority selection.} 

This third event happens with probability p. 

Thus the process can be coded by sequences of type {-Knew? -Roid? -Rcomp}^- 

Let Xt be the t-th element of this sequence. Then the renewal at every i?oid event can 
be formalized as follows: Using the fact that the sequence {Xt}^i is independent of the 
priorities, for every x G M 

F{Ot+i <x\Xt = Rom) = <x\Xt = i?oid) = F{Nt < x) = 

where recall that Ot and Nt stands for the priority of the old and the new task at time 
step t, respectively. 

Since the priority distribution of the old task is renewed at every event -Roid, in order 
to determine the priority distribution of the old task in the system for a big enough T , 
we have to count how many -Rcomp-s we had after the last i^oid event. 
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It is easy to see that the distribution of the arrival time of the last i?oid event is just 
T — GEO (^^) + 1- Similarly, if we only want to count the number of -Rcomp-s after 
the last -Roid; we have to re-normalize the probabilities to exclude cases of -Rnewj so it's 
distribution is 

GEO I I - 1 = GEO ( ^ ] - 1. 



Thus, the old task in the system at time T had to compete with the other tasks at 
each event -Rcomp and had to "win", i.e its priority had to be less than all of them, in 
order to stay in the system up to time T. So, its priority had to "defeat" GEO - 1 

many other tasks to stay in the list so far. Hence the distribution we were looking for 
is the minimum of X ~ GEO (^j^^ many independent random variables ({l^jj^i), of 
distribution R{x). Let us compute it's CDF using the law of total probability: 

oo 

F{x) = P( min Yi < x) = ^P(3i < X : Yi < x \ X = k)F{X = k) 
*~ ■■' fc=i 

1 + P A 1 



1 + 

in agreement with [I I], since the summation converges if p < 1. A comment is that we 
only used the extension to doubly infinite stationary sequence since the geometric variable 
in this proof is unbounded. 

3. Records and expected waiting time 

3.1. Ergodicity for p < 1. In this section we verify Conjecture 1.1 which gives us the 
expected waiting time in the Barabasi case. Let us first set p < 1 and L > 2. We use 
the method suggested by Vazquez [l i], and work out the details. Using the arguments 
introduced by Vazquez let us denote by Tj the waiting time of the task executed at step 
i, and by r/^ (/ = 1, . . . , L — 1) the resident times of the tasks that are still in the buffer 
at step t. Then we summarize the main point in the following theorem: 

Theorem 3.1. [Ergodic Theorem for L > 2,p < 1] The priority queuing system in the 
Barabasi case is ergodic for any L > 2 and p < 1 and the following limit holds: 

1 * 

(3.1) IE/^(t) = lim - \ Ti = L a.s. and in Li. 

Proof. Since the total buffer time until time t is Lt and in each time step exactly one task 
leaves the system, we have that 

t L-l 

^ + ^ r/t = Lt. 

i=l 1=1 
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Hence 



1 * 1^"^ 
(3.2) lim-5^r, = L- lim-^r/, 

i=l 1=1 

if these limits exist. Similarly to the L = 2 case in Section 2, the process renews whenever 
there are L — 1 successive Rom events, where i?oid means that the oldest task is selected 
for execution because of the random selection. For every t the probability that there is a 
renewal section of the process starting at the step t is {^-j^)^ ^ ■ In addition, two renewal 
sections starting at time ti and ^2 are independent if \ti ~ > L — 1. So the system 
renews infinitely often a.s.. (This implies also the existence of a stationary distribution 
for the L — 1 not newly added task's priority.) Analogous to the argument in Section 2, 
each remaining time r/^ can be stochastically dominated by a geometric random variable 
with success probability {^-j^)^"^- Thus, ^ — t- as t — t- oo holds for any I almost surely 
(a.s.) and in Li, and so we obtain 



(3.3) lim -} t[, = a.s. and in Li. 

1=1 

There remains the proof of ergodicity. Let us denote by M_^ = {M^, . . . ,M^~^) the 
L — 1 dimensional vector of the oldest elements in the buffer at time t, such that is 
the priority of the oldest, is the priority of the z-th oldest task in the system. Then 
clearly the sequence {M_^}'^-^ is Markovian with state space [0, l]^"-*^. Moreover, it is 
clearly aperiodic, and since the system renews infinitely often a.s., it is irreducible. Thus 
ergodicity and equations (3.2) and (3.3) immediately imply (3.1) as it was suggested by 
Vazquez's heuristic arguments [l I]. □ 

3.2. Records. In this section we investigate the case p = 1, L = 2 , i.e. the case when 
the selection protocol has no randomness and always picks the task in the system with 
higher priority. We show that the system achieves no stationary distribution in this case. 

It is easy to see that when a task arrives whose priority is less than all the priorities 
before, then the old task in the system will be chosen for execution and this task will 
remain in the buffer. This implies that Ri{x, t) (the distribution of the old task's priority 
at time t) is the distribution of the minimum, i.e. the lower record value of the first t + 2 
tasks' priority. That is, we need to consider the minimum of t + 2 i.i.d. random variables 
with distribution R{x). 

On the other hand, a stationary distribution Rl{x) should satisfy equation (1.4), which 
has only a degenerate solution Ri{x) = a.s. if p = 1. Thus, the stationary distribution 
Rl{t) will never be achieved, since 

P(Vt Ri{x,t) 7^ 0) = 1. 

Hence, the properties of the waiting time r are related to the process of records. 
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In order to introduce some classical theorems of record theory (for an introduction see 
[1, p. 22-28]), we will need some definitions. 

Definition 3.2. Let {Xt}^^ be i.i.d. random variables with absolutely continuous CDF 
F . Then we define the fallowings: 

(1) The sequence {T^}^^^ of lower record times is defined by 

Ti = 1, and for k > 1 

Tk = min{t | Xt < Xt^_^}. 

(2) The sequence {Afej^j of inter-record times is defined by 

= Tk — Tk-i- 

(3) The sequence {It}tLi of record indicators is defined by 

/i = 1 , and for t > 1 

/t = l{Xt<min{Xi,...,Xi_i}}. 

(4) The sequence {xk}'k=i of lower record values is defined by 

Xk = Xt^. 

If we do not care about the record values, without loss of generality one can assume that 
the distribution of the Xj-s is f/[0, 1]. Thus one can easily compute the joint distribution 
of the If-s. Let n, 1 = ti < ^2 < • ■ • < be arbitrary positive integers. Then 

P(/,, = l,...,/i„ = l) 

= / P(/ti = 1, . . . , /t„ = 1 I Xt,=^^,. . . , Xt„=^Jdxi ...dxn 

J0<Xl<--<Xn<l 

— / ™<2-l™t3-*2-l ^t„-t„-l-lj J 

— I 1 2 ' * * n Ujuj\ . . . Ujdj^ 

Jo<xi<--<x„<l 
_ 1 1 1 

tl t2 tn 

Since the integers n,l = ti < t2 < ■ ■ ■ < tn were arbitrary, and the It-s are discrete random 
variables, it is evident that these random variables are independent, and It ~ BER 
is a Bernoulli random variable. 

Now consider the ratio The distribution of ^ conditioned on the value of Tk-i can 
be determined by Tata's reasoning [13]: 



P ( — > s 



Tfe_i = t] = F{ ^\ > X 

P ( Afe > -^t 

1 — X 



Tk-i — t 



Tk-i — t 



P(/i+i = 0,...,/LV(i-.)j =0) 
t 
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Since — )■ oo a.s. as A; — t- oo, hence 

lim P ( ^ > x] = 1 — X, 

SO thus jr- and are asymptotically U[0, 1]. 

Now we turn back to the Barabasi model in the case p = 1 and L = 2. Equation (3.2) 
holds also in this situation: 



1 * '^[t 
lim - > Tj = 2 — lim — - 

t— >co t t-^oo t 

i=l 

if these limits exist. But now, with the help of record theory one can prove that the limit 
on the right hand side does not exist. On the one hand if we choose t = + 1 then we 
obtain: 

liminf — - = a.s., 

t—^oo t 

since the sequence of records is infinite a.s. (This is ensured by the continuity of the CDF 
F), so thus = will hold for infinitely many t-s. On the other hand the limsup is 
taken if t = (we look the subsequence of records): 

hm sup — - = lim = hm 

t-s-oo t fc-s>oo Tk k-^oo Tk 

= 1 _ lim %i = 1 - Z, 

A;— >oo Tk 

where Tk and are the fc-th record and inter-record times, and Z is a f/[0, 1] random 
variable. 

Thus, Vazquez's heuristics does not work in this case. Note that the dynamics is not 
even ergodic, since the irreducibility does not hold, as the priority of the old task decreases 
monotonically. 

Further, we can state precise asymptotic results about the waiting time of the records 
in the queue. It was proven by Holmes and Strawderman [8] that the Strong Law of 
Large Numbers (SLLN), by Neuts [ ] that the Central Limit Theorem (CLT) and by 
Strawderman and Holmes [ ] that the Law of Iterated Logarithm (LIL) hold for the 
series of In Tk and In A^. as well: 

Theorem 3.3. Let Zk he either Ak orTk- Then 
(1) [SLLN][ ] 

InZfc a.s. ^ , 

— )■ 1 as k ^ oo. 



k 

(2) [CLT][ ] 

\n.zi. — k 
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(3) [LIL][12] 

T { In Zk-k \ 

limsup< = ( = ^ d-S-, 

fc-5-oo [_ \/2k\n\nk ) 

liminf | ^ = \ = —1 a.s. . 
k^oo V2A;lnln A; J 

Thus we see that the k-th record spends roughly e'^ time in the buffer. 

4. A NATURAL GENERALIZATION OF THE MODEL 

In this section we consider a common generalization of Barabasi's model and Example 
1.2 which were not analyzed in the literature before. Also in this new model there is a 
list with L tasks, the new tasks has i.i.d. non-negative priorities, and in each discrete 
time step one task is selected to execution and a new task arrives to it's place. The only 
change is that the dynamics of the system is now given by a general selection protocol. 
We will analyze this new model for the case L = 2, so we define the selection protocol to 
this case. In every time step we select a task according to the priorities in the following 
way: 

(4.1) v{x,y) := P(new is chosen | new has priority x, old has priority y), 

where 

Assumption 4.1. We assume that "^y v{.,y) is increasing and^x f (x, .) is decreasing. 
Furthermore let us suppose that 

Wx Wy v{x, y) < ci < 1 

with some constant Ci. 

Let Ot be the priorities of the old task at time t. Clearly, the sequence {Ot}'^^ is an 
aperiodic Markov-chain on the state space (0, 1). 

Lemma 4.2. Under Assumption 4-1, the Markov-chain {Ot}'^i defined above is positive 
recurrent. 

First note that positive recurrence of the Markov chain implies the existence of a station- 
ary distribution as well. We remind the reader that q{s) and qi{s) was defined in (1.2) and 
(1.3), and r{x) is the density function of the new task. Since 3ci Vx v{x,y) < Ci < 1, 
hence Vs q{s) = v{y, s)dR{y) < Ci < 1 and Vs ^1(5) > 1 — Ci > 0. 

Proof. For any x G supp{r) let us denote by T^^y^e the hitting time from x to {y — e,y + e). 
Thus the distribution of T^^y^e is stochastically dominated by the geometric distribution 
with success probability (1 — ci){R{y + e) — R{y — e)). Therefore 

^{Tx,y,e) < JZ ^7^7 ^7 ^ < 00. 

(1 - Ci)iRiy + e) - Riy - e)) 

□ 
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Since the Markov chain is ergodic, this fact imphes that Theorem 3.1 is vahd in this 
case, and the expected waiting time of a task equals 2. 

4.1. Distribution of the priorities. In the following we write the density function ri(x) 
of the old task's priority at stationarity as a fix point of a Hilbert-Schmidt operator. Then 
using this fix point equation we show an opportunity of approximating ri{x), and then 
we approximate ri{x) of Example 1.2. 

Let us consider the general case. We know from Lemma 4.2 that stationarity distribu- 
tion exists for the chain. The stationary equation (1.4) is 

Ri{x) = / ri{s)q{s)ds + / r{s)qi{s)ds. 
Jo Jo 

After substituting (1.2) and (1.3) we obtain 

PX pi 

Ri{x) = / {ri{s)r{y)v{y,s)+r{s)ri{y){l-v{s,y)))dyds 
Jo Jo 

Differentiation with respect to variable x yields the following equation for the stationary 
density function ri(x) of the old task's priority: 

ri{x) = / {ri{x)r{y)v{y,x) +r{x)ri{y){l -v{x,y)))dy 

Jy=0 

With the constant Ci from the definition of w (4.1) we can write: 

ri{x)(^l-j r{y)v{y,x)dy^ = r{x) i^j ri{y){l - v{x,y))dy - + Cir{x). 
It is easy to see that the coefficient of ri(x) on the left side cannot be zero, hence 



, r{x){l-v{x,y) -ci , 

n{x) = I -I — —n{y)dy 

1 - J_,^Qr{z)v[z,x)dz 



y= 



(4.2) 

^ Cir[x) 



1 - J^^Qr{z)v{z,x)dz 
To shorten the notation we introduce 

a{x,y) := 1 - v{x,y) - ci, 
r{x) 



9[x) 



1 — r, r(z)v(z, x)dz 
(4.3) •^^=° ^ ^ ^ ' ^ 

K{x,y) := a{x,y)g{x), and 



1 — /^=o x)dz 
Moreover we introduce the integral operator A associated to kernel function K, i.e. 



Aip{x) := j K{x,y)ip{y)dy. 

Jy=0 
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Since K G ^^[0, 1] x [0, 1], hence A maps L^fO, 1] into itself (see [ , Theorem 9.2.1.]). With 
this notation, equation (4.2) can be written as 

(/-A)ri(x) = /(x). 

Hence, assuming that the operator I — A is invertible for certain choice of the underlying 
distribution R and the selection protocol v, we would like to determine the function ri(x) 
which satisfies 

(4.4) n{x) = {I-A)-'f{x) 

If X]^o^"/('^) converges we could write 

oo 

nix) = J2Ay{x). 

n=0 

It turns out that A^ is also an integral operator with some kernel function Kn, which is 

"1 pi 

Kn{x,y) = I ... I K{x,ui)K{ui,U2) . . . K{un-i,y)dui . . .dUn-l- 




Note that we can benefit from the fact that f{x) = Cig{x) (see equation (4.3)), so we get 



A"/(x) = A"^(x)ci = ci / K^{x,y)g{y)dy 



ci / ... / a{x,ui)g{x)...a{un-i,t)g{un^i)g{y)dui...dun-idy 



JO 



cig{x) I ...I a{x,ui)g{ui) . . .a{un-i,Un)g{un)dui. . .dun- 



Jo 



K{x,u\) 



Thus we see that it is useful for us to introduce another kernel function given by 

(4.5) K{x,y) := a{x,y)g{y), 
and the corresponding integral operator 

(4.6) Ah{x) := [ K{x,y)h{y)dy, 



which maps L^[0, 1] into itself. With this notations, the terms A"/(x) can be obtained as 
cig{x)Hn{x), where 

Hn{x) = I ... a{x,ui)g{ui) . . .a{un-i,Un)g{un)dui . . .dUn 
Jo Jo 

... / K{x,Ui) . . . K{Un-l,Un)dUi . . . dUn = A^'i{x). 
10 JO 

So we see that (4.4) is equivalent to 

(4.7) ri(x) = ci g{x){l + H^{x) + H^ix) + ...), 
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which yields: 

ri(x) ^ ci g{x){l + Hi{x) + H2{x) + ■■■ + Hn{x)). 

4.2. Example. In the following we consider Example 1.2. Recall that in this case, priority 
service happens with probability and in this case a task is chosen proportional to its 
priority, i.e. v{x,y) = px/{x + ?/) + (! — p)/2. In the sequel, using the software Wolfram 
Mathematica, we determine ri{x) numerically for this particular choice with some fixed 
parameters. We compute the terms in (4.7) and give estimates on Hn{x). From these 
we can derive estimates on parameters {p, c) where the series in (4.7) converges. Before 
doing these, we argue shortly that shifting the distribution of the new task's priority up 
by c has a real relevance in the literature. To see this, note that if X, F ~ U[c, 1], and 
V,W ^ t/[0, 1], then 

X V{l-c) + c V + 5 

X + Y V{l-c) + c + W{l-c) + c V + W + 26' 

where 6 = > 0. This protocol is similar to the growth rule of the preferential 
attachment model, see [3, 5, 4]. 

Recall the definition of the kernel function K and the corresponding operator A. We 
have seen that Hn{x) = 74"l(a;). So in order to prove convergence of (4.7) it is enough to 
give an upper bound on the norm of if„(x). This can be achieved by estimating the 
norm of A or the Hilbert-Schmidt norm of K: 



\HJx)h < \\AK < 



where ||-ft'|||^5 — £ £ K'ix,y)dydx. If \\K\\ 

Hs ^ ^ then the sum in equation (4.7) 
converges, so it remains to find the parameter region (p, c) for which this holds. Using 
the definitions given in Example 1.2 and in equations (4.3) and (4.5) we get that 



l^ll^^ = / / I ^tt:; i±p _ h In 1+^) —c 1 ^^^^ 

c+x / . 




2 



^In^) 1 

— c c+x J 



2 

dy. 



The expression increases in p since in the denominator decreases. This implies that 
once we have found a pair (p, c) for which < 1 then this inequality also holds for 

all pairs (g, c) where Q < q < p. The region of pair of parameters for which ||A|||/5' < 1 
can be determined numerically, and it is represented on Fig. 1. In addition. Fig. 2 and 3 
illustrates the density ri(x) and the distribution of r on a log-log plot for three different 
choice of parameters, respectively. 
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Figure 2. The DF of the old task's priority in the stationary case with 
c = 0.2 and different values of p. Blue: p = 0.7. Red: p = 0.8. Mustard: 
p = 0.9. 
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Figure 3. The distribution of r, the waiting time of a new task in the 
stationary case, with c = 0.2 and different values of p drawn on a log-log 
plot. Blue: p = 0.7. Red: p = 0.8. Mustard: p = 0.9. 

Now we give estimates on the distribution of r when k > 1, which is given in equation 
(1.5). In this example the functions q{x) and ^1(0;) (defined in (1.2) and (1.3)) can be 
written in the following form: 

p y 1—p l+p px 1 + X 
q{x) = / ay H — = — m 



1 — c X + y 2 2 1 — c c + X 

y 1 — p 

Qiix) = p / — ^ — dRi{y) ' 



X + y 2 
It is easy to see that q'{x) 7^ 0, so we can write 



1 

fc-2. 



{T = k)= I qi{x){l - q{x))q{xY-'dR{x) 

-{-q'{x))q{xY~'^dx. 



^ 1 qi{x){l - q{x)) ^ „//^^^„/^^fc-2, 



1 — c —q'[x) 



■.=L{x) 

With the new notation L[x) := ^^^^^l^^^f^' claim that L(x) can be bounded 

from below and from above by some constants m, M only depending only on {p, c): 



(4.8) 
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Since 

1 — p 

— \-pc< gi(l) < qi{x) < 1, 

1 -p 1 

— h pcln — < 1 - g(c) < 1 - q{x) < 1, 



Thus if we set 



1 - C , 1 . 1 . 1 {P,c)^(l,0) , ^ 1 
— Ti- < < < > m 2 H — . 



-..:=-^#^^f#, and Mp,c:^ ^ 



(1-cMl) ' • (l-c)g'(l)- 

then the bound (4.8) follows. With these bounds, one can easily estimate P(r = A;): 



1 1 



{T = k)< M{p, c) I q'{x)q{xy-^dR{x) = M{p, c)- -{q^-\c) - ^'^"^(l)) 

K, 1 



and similarly 
In the limit (p, c) — )■ (1, 0) we get 



T = k)> m{p,c)^{q'~\c) - q'-\l)). 




(p,c)^(l,0) 

Since g(c) > g(l), hence 

(4.9) P(r = k)^ const ■ \q^{c) = const ■ \e~^'^\ 

where 

= 

ln(g(c)) 

This means that the probability that the waiting time is of length k behaves approximately 
as 4 until k < kQ, i.e. it obeys a power law-like behavior. Furthermore, if c 7^ and 
{p,c) — )■ (1,0), then ko 00, i.e. the exponential cutoff in equation (4.9) shifts to 
infinity, and thus for {p, c) values close to (1,0) the distribution of r will be close to a 
power-law distribution. 

5. Summary 

In this paper we investigated and generalized the priority queueing model of Barabasi. 
We showed that in the original model of Barabasi, the system is ergodic and irreducible as 
a Markov chain once the priority selection probability is separated away from 1. Further, 
we gave a more probabilistic approach to compute the distribution of the priority of 
the old task in the system if the buffer length equals 2, and determined the average 
waiting time of a task in the system for arbitrary buffer length. We investigated the 
case p = 1 separately and found that the model is equivalent to the processes of records. 
Next, we generalized the priority queueing model with an arbitrary selection protocol 
depending on the priorities in the system and some extra randomness. We found that 



18 JULIA KOMJATHY, KAROLY SIMON, AND LAJOS VAGO 

the system is ergodic and irreducible once the selection protocol is separated away from 
1. Further, we gave a description of the density of priority of the old task in terms of 
Hilbert-Schmidt operators. We investigated a special example when the priority selection 
is done proportionally to the priorities of the tasks in more detail. Namely, using bounds 
and approximation methods based on the operator approach, we gave a region of pairs 
(p, c) where the Hilbert-Schmidt description of the density function is converging and 
determined the waiting time distribution in this case. 
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